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Galactic HI on the 50- AU scale in the direction of three 
extra- galactic sources observed with MERLIN 
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ABSTRACT 

We present MERLIN observations of Galactic 21-cni Hi absorption at an angular 
resolution of ~ 0.1 — 0.2 arcsec and a velocity resolution of 0.5 km s^^, in the direc- 
tion of three moderately low latitude (-8° < b < -12°) extragalactic radio sources, 
3C 111, 3C 123 and 3C 161, all of which are heavily reddened. Hi absorption is ob- 
served against resolved background emission sources up to ~ 2 arcsec in extent and 
we distinguish details of the opacity distribution within 1-1.5 arcsec regions towards 
3C 123 and 3C 161. This study is the second MERLIN investigation of small scale 
structure in interstellar H I (earlier work probed Galactic H I in the directions of the 
compact sources 3C 138 and 3C 147). The 0.1-arcsec scale is intermediate between 
H I absorption studies made with other fixed clement interferometers with resolution 
of 1 to 10 arcsec and VLBI studies with resolutions of 10-20 milli-arcsec. At a scale 
of 1 arcsec (about 500 AU), prominent changes in Galactic Hi opacity in excess of 
1-1.5 are determined in the direction of 3C 161 with a signal-to-noise ratio of at least 
lOtT. Possible fluctuations in the Hi opacity at the level of about 1 are detected at the 
2.5 — 3cr level in the direction of 3C 123. 

Key words: radio lines: ISM - radio continuum: galaxies - ISM: kinematics and 
dynamics - ISM: structure - Galaxy: kinematics and dynamics 



1 INTRODUCTION 

In order to understand the dynamics of the interstellar 
medium (ISM), a thorough understanding of the physical 
conditions of the ISM gas is required over a large range of 
physical scales. In particular, the magnitude of changes in 
density of the neutral atomic component (H l) over a wide 
range of physical scales from tens of AU to parsecs plays a 
major role in determining the evolution of the dense compo- 
nent of the ISM. A recent conference in Socorro, New Mexico 
(USA) in May 2006 was devoted to this broad topic of small 
scale structures in various ionized and neutral (atomic and 
molecular) phases of the ISM (the SIN S conference: 'Small 
Scale Ionized and Neutral Structures'; iHaverkorn fc GossI 
120071 ). A major topic of discussion at this conference was 
the nature of Hi structures with scales from tens to thou- 
sands of AU. To date, most of the effort in studying the small 
scale structure in H I has consisted of V LBI and pulsar de- 
terminations at the 10 to 100 AU scale (Faison et al. 199S; 
Faison fc Goss 2001; B rogan et a l. 2005; Brogan et al. 2007; 



Weisberg fc Stanimirovicll2007l : IStanimirovic et aLlbOOd l . 

A major gap between the 100- AU and O.l -pc (2 x 10"^ 
AU) scales has remained largely unexplored; iDavis et al.l 
j 19961 ') used the MERLIN array and the European VLBI 



Network (EVN) to observed Galactic H l absorption towards 
the extragalactic sources 3C 138 (Galactic coordinates I = 
187.4° and b = -11.3°) and 3C 147 (Galactic coordinates I 
= 161.7° and b = 10.3°) over angular scales of 50 to 1000 
milli-arcsec (mas). Prominent changes in Hi opacity at the 
level of 0.15 were observed in the MERLIN 0.1-arcsec reso- 
lution data and at the level of 0.5 in the 50-mas EVN data. 
The changes were observed over linear scales of < 500 AU 
arising from nearby Hi gas (< 500 pc). 

A number of discussions of the observational charac- 
teristics of the small scale neutral gas in addition to H I 
were presented at the SINS conference (JLauroesch 1 120071 ) . 
[Andrews et al.l (J200ll ) have reviewed the nature of both 
spatial and temporal v ariations in optical IS lines, while 
Ijenkins fc Trippj (120071 ) have reviewed C I UV observations 
suggesting overpressure in some portions of the ISM. If the 
observations of opacity variations are interpreted in terms 
of density fluctuations, a major problem of excess density 
and p ressure is implied. To ove rcome this problem iHeilesI 



11993) (see also lHeilesllioOTi and lHeiles fc Stinebring||2007l ) 
has suggested several geometrical explanations that would 
reduce the ov erpressure; exampl es includ e sheet or filament 
morphologies. lOeshpandel l|2000h (also see lOeshpande 1120071 ) 
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has, on the other hand, proposed that the structures on 
small scales of tens of AU are only a natural extension of 
irregularities in the distribution of H I scales at much larger 
di mensions, based on a po wer law ('red spectrum') derived 
bv lDeshpande et alj (|2000l ) using VLA data in the direction 
of the Galactic source Gas A. This derivation consisted of 
a power law over scales of several parsec to less than 0.1 
pc. The pro s and cons of these contra sting proposals are 
discussed bv iHeiles fc Stinebriiigl (|2007l ). 

In order to fill in the gap between the 0.1-pc and 10- AU 
scales, additional observations of the 100- AU scale are re- 
quired. The MERLIN array at 21 cm is the ideal instrument 
for this endeavour; we undertook H I absorption observations 
in the direction of strong extragalactic sources with total 
angular sizes of order 1 arcsec. Three sources with Galactic 
latitudes in the range -8° to -12° were observed in 2001 and 
2003. These observations repr esent the second in stalment of 
Galactic H I studies following iDavis et al.l l|l996r ) . 

The choice of sources was guided by the pioneering (and 
under- reco gnized) ground-br eakin g Galactic Hi absorptio n 
studies bv iLiszt et all l| 19821 ) and iGreisen fc Liszt I (|l986h . 
Observations were made in the early 1980s, using the newly 
completed VLA with an angular resolution of 3-4 arcsec and 
a velocity resolution of 1.3-1.6 km s~^ , in the directions 
of 3C 111, 3G 123, 3C 161 and 3C 348. The new MER- 
LIN observations were made with an angular resolution of 
~ 0.1 — 0.2 arcsec and a velocity resolution of 0.49 km s~^, in 
the direction of all of these sources except for 3G 348. Pre- 
limina ry MERLIN results were presented by iFaison et al.l 
i2003). 

We outline the data acquisition process in Section [5] 
and describe and interpret the Hi absorption towards each 
source in Section |3] The results and implications for the 
structure and kinematics of Galactic H I are summarised in 
Section |4] and we add an Appendix |X] on hitherto unpub- 
lished data for 3G 161. 



compact phase reference sources listed in Table [T] We used 
a 0.5-MHz bandwidth divided into 256 channels for the H I 
targets and the bandpass calibration source, giving a veloc- 
ity channel separation of 0.41 km s~^. The target data were 
recorded adjusted to a fixed velocity with respect to the local 
standard of rest (VLsr), given in Table [T] Other calibration 
sources were observed using the standard wide-band contin- 
uum configuration for maximum sensitivity. The absolute 
flux density scale for all sources was derived from 3C 286. 
We used 3G 84 to calibrate the bandpass, after eliminat- 
ing its weak Galactic absorption. We applied these solutions 
and the phase and amplitude calibrations derived from the 
phase reference sources to the relevant targets. More de- 
tails of MERLIN spe ctral line data reduction are given in 
[Diamond et all l|2003D . 

After calibration, we inspected the spectrum of each 
target and identified the line-free channels. These were sub- 
tracted from the spectra to enable us to make H I absorption 
line cubes and we also imaged the continuum. The final op- 
tical depth cubes were calculated from the continuum image 
and the line cube. We used cutoffs of a few per cent of the 
peak continuum intensities (see Table [T] and captions to the 
relevant figures). All positions are given in J 2000 (ICRS). 
We can only measure Hi absorption on the same scale as 
the continuum structures, since H I emission is undetectable 
on the (sub-)arcsec MERLIN scales. The maximum scale of 
emission which can be imaged reliably from our observa- 
tions is 2-3 arcsec at 1.4 GHz and 0.7-1 arcsec at 5 GHz. 
Detectable emission on slightly larger scales would cause a 
distinctive 'bowl' artefact due to missing flux density on 
most baselines. No signs of this are detected in any of the 
MERLIN images, implying that there is minimal missing 
flux density on the scales from 0.1 to 3 arcsec imaged here. 
In addition we compare the current MERLIN results with 
existing VLA data to conflrm the minimal effect of missing 
short spacings. Detailed comments are presented for each 
source. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were carried out in a two-stage process. 
The structure of our targets had never been investigated 
previously at 0.1-arcsec resol ution (with the excep tion of 
the 3G123 hotspots at 5 GHz. iHardcastle et al.r l997) so we 
started by observing the radio continuum from each source, 
at a wavelength of either 6 or 21 cm, in order to determine 
the size of the field required for the final H I line and contin- 
uum images. The continuum observations were made using 
the MERLIN array including the antennas at Defford, Cam- 
bridge, Knockin, Darnhall, Tabley, and the MK2 telescope 
at Jodrell Bank, giving a usable bandwidth of 13-15 MHz. 
The observing parameters are given in in Table [T] Stan- 
dard observing and data reduc tion procedures were used, 
see the MERLIN User Guide (JDiamond et al.l 120031 ). The 
resulting continuum images are in good agreement with the 
somewhat lower-sensitivity images made during the H I ob- 
servations (see below). The continuum emission from 3C161 
is discussed further in Appendix El 

The observing parameters for the Hi observations are 
given in Table [T] The antennas listed above were used with 
the addition of the Lovell Telescope. Observations of the 
spectral targets were interleaved with observations of nearby 



3 RESULTS 



3.1 3C 111 



The large low latitude source 3C 111 {I = 161?68, h = -8?82) 
has a complex continuum structure. This FR (Fanaroff- 
Riley) II source is a nearby broad line radio galaxy at a red- 
shift of 0.049. VLA continuum studies wer e carri ed out by 
iLindfield fc PerlevI (| 19841 ) : lGreisen fc Lisztl (|l986l ) have also 
summarized VLA H I observations of this source, identifying 
a core continuum component (their '2' ) and two prominent 
lobes with hotspots, NE (their '1', displacement 2 arcmin) 
and SW (their '3' and '4', displaced by 80 arcsec). Using 
the MERLIN array, the NE component is detected with an 
extent of ~ 2.5 arcsec and a continuum peak of only 46 
mjy beam~^ (not listed in Table [T|; the mean Hi spectrum 
was obtained after smoothing to a 0.3-arcsec resolution. The 
core of the galaxy has a peak of 1.27 Jy beam~^ . The two H l 
spectra for components 1 and 2 are superimposed in Fig. [l] 
The SW component is heavily resolved by MERLIN; 
the emission detected over ~ 2 arcsec has a peak of only 8 
mJy beam~^ . At this level, detection of H I absorption is not 
possible even after some smoothing. We also note that the 
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highest apparent H I opacities for the NE component (at -5 
and km s~^) are unrehable with uncertainties in excess of 
1, due to the partial saturation of the H I line and the weak- 
ness of the resolved continuum feature. The faintness of the 
lobes of 3C111 implies that it is impossible to determine de- 
tailed H I opacity images. The spectra obtained for the core 
source and for the average of the NE lobe using MERLIN 
are in excellent agreement with earlier VL A data obtained 
with a resolution of 3.5 arc sec obtained bv lGreisen fc Liszt I 
l|l986r ). We confirm the increase in the H I opacity for the -19 
km s~^ feature, from r Hi ~ 0.35 towards the core (compo- 
nent 2), compared with r Hi < 0.1 towards the NE hotspot 
(component 1). A relatively higher opacity of r Hi ~ 0.6 is 
observed towards the NE hotspot at -11 km s~^, compared 
wit h T Hi ~ 0.3 toward s the c ore. 

iMoore fc Marscherl ([1995|) made observations towards 
3C 111 using the VLA at 16-14-arcsec resolution at the re- 
spective Hi and OH line frequencies of 1.4 and 1.6 GHz 
as well as multi-epoch observations of H2CO at 4.8 GHz 
in various configuration s. Their images are als o consistent 
with the Hi results from lGreisen fc Liszt I (|l986l ). containing 
substantial absorption line variations in Hi and OH which 
suggest with structure on scales of 0.14 to 0.34 pc at their 
adopted distance of 350 pc. The time variability in H2CO 
absorp tion observed towards the core by Moore fc Marscherl 
(|l995l) and iMarscher. Moore fc Banial |l993h over a 3.4 vr 
period provides additional evidence for 10-AU concentra- 
tions with number densities > lO'' cm^^. These apparent 
fluctuations could arise from a combination of the elliptical 
apparent motion due to paralla x, linear solar motion an d 
the proper motion of the cloud (jMoore fc Marscherl 1 1995h . 
giving a transverse drift of about 4 AU per year (at the 
adopted distance of 350 pc) of the position of the line of 
sight through each cloud towards the background source. 



tions. The narrow line at -20.1 km s"^ (FWHM 2.9±0.2 km 
s"^) and the +4.2 km s'^ feature (FWHM 6.0±0.5 km s~^) 
are prominent; the latter has been associated with a CO line 
indicative of a possible association with a molecular cloud 
(JLiszt et al.lll982l : lLiszt fc Burtonlll979l '). The mean opacity 
of the -20.1 km s~^ line in the M ERLIN data is som ewhat 
larger than the value reported bv lLiszt et al.l (|l982l ) (their 
position 'c'); the properties of the -(-4.2 km s~^ feature agree 
closely. In the absence of a direct measurement, a possible 
distance d is given by 



d =< 1^1 > /sinfe. 



(1) 



where z = 107 p c is the mean value of the Galactic height 
(|Crovisierlll978l ). This suggests that the moderate- velocity 
H I observed tow ards 3C123 is at d ~ 530 pc (as adopted by 
iLiszt et al.ll 19821 ). where 0.1 arcsec corresponds to transverse 
sizes « 50 AU. The absorption at -20.1 km s~^ probably 
arises from more distant gas but a more precise estimate is 
not possible in this direction. 

In Figs.|4]and[5]we show the major new results based on 
the MERLIN Hi data for 3C 123. Figure g] shows the optical 
depth image for the strongest line at 4.2 km s~^. At this 
velocity, variations in the Galactic Hi opacity are observed 
with a statistical significance of only 2 — 3a". The signal-to- 
noise ratio is limited due to the low surface brightness of 
this resolved source, using a beam size of ~ 0.14 arcsec. The 
angular scales of the variations vary from 0.2 to 0.6 arcsec 
(100 to 300 AU). In the N-S slice a variation of Ri 1 ± 0.4 in 
opacity is observed. In the E-W slice a variation in opacity of 
0.7 ±0.25 is observed on a scale of 0.2 arcsec and a variation 
of 1.5±0.5 over 0.3 arcsec. Figure [S] shows a comparable plot 
of the variation of Hi opacity at -20.1 km s~^. The weak 
absorption is essentially constant over the scale of ~ 2 arcsec 
(1000 AU) at a level of r (Hi) ~ 0.1. 



3.2 3C 123 

The continuum image of 3C 123 {I = 170?58, b = -11?66) is 
shown in Fig. [5] Previous high- resolution VLA images of this 
FR II radio galaxy have been presented by 1 Cox et al.| (|l99l|) 
(images made by R. La ing at 1.7, 4.9 and 14.9 GHz) and by 
iHardcastle et al.l l|l997i ). 3C 123 has been identified with a 
19 mag galaxy with strong Gal actic reddening at a redshift 
of 0.22. IHardcastl e et al.l l|l997h show a 5 GHz MERLIN im- 
age (beam 0.05 arcsec); both components shown in Fig. [5] 
(the E and W hotspots) have been imaged by these authors. 
The weak core (~ 50 mjy) is also detected in the current im- 
age at right ascension 04''37'"04'!38, declination 29°40'13'.'8 
(J2000). The MERLIN continuum image of 3C123 contains 
all the fiux density that was detected by the VLA at 1.4 GHz 
with 1.2 arcsec resolution (R. Laing in Atlas of DRAGNS 
lhttp://www .jb.man.ac.uk/atlas/ ). The total fiux density of 
13.5 Jy above 3a in the MERLIN image of the E hotspot 
(Fig. 12)) in an area of about 10 arcsec"^ is in good agreement 
with the VLA image over the same region. The central core 
radio source is unresolved by MERLIN and has a fiux den- 
sity of 48 mJy, compared to 51 mJy in the VLA image. 

The mean Hi optical depth spectrum toward the SE 
hotspot (over an area of 0'.'25 x 0'.'25) is shown in F i gure |3l 
This spectrum is similar to the results of lLiszt et al.l (|l982h . 
Their anomalous velocity feature at -73 km s~^ is outside 
the velocity range (± ~50 km s~^) of the current observa- 



3.3 3C 161 

This low-latitude heavily reddened field {I = 215? 44, b = 
-8? 07) remains an empty field with no certain optical iden- 
tification. 3C 161 is comp act (|Ulvestad et al.l|l981 ; PerlevI 
ll982l : |Pearson et al.lll985h with observed radio structure on 
the arcsec scale. The later authors present a VLA image 
at 5 GHz (be am 0.4 arcsec) show i ng an obvious extension 
to the West. 'Edwards fc Tingavl l|2004h have pointed out 
that 3C 161 is a CSS (compact steep spectrum source) with 
a spectral index of ~ —1.2. The source has also been dis- 
cussed by iKing et al.l (|l996l ) . The continuum im age made 
during the VLA Hi observations of iGreisen fc L iszt (1986|) 
(beam 4.2 arcsec) remains a remarkable contribution in pro- 
viding the overall morphology of the fainter components of 
the source over a range of ~ 1 ar cmin (see Appendixf^)) . The 
source 'B' 2 arcsec to the NE (lUlvestad et al.lll98ll ') is not 
detected bv lGreisen fc Liszt 1 1| 19861 ). However, the strongest 
compo nent of 3C 161 has not been detected by VLBI tech- 
niques; |Fomalont_et_alJ (|2000r ) did not detect 3C 161 during 
the VSOP prelaunch VLBA 5 GHz survey. Faison and Goss 
have also observed 3C 161 with the VLBA in 1998 (unpub- 
lished) with a faint detection of emission on an angular scale 
of tens of mas. Thus 3C 161 contains no prominent milli- 
arcsec scale structure and is an ideal candidate for MERLIN 
0.1 arcsec Hi imaging. 3C 161 is the brightest of our three 
sources (the continuum fiux density exceeds 0.1 Jy beam~^ 
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over almost 1 arcsec^ ) , giving very good sensitivity to H I 
absorption on scales of 100 mas. 

The continuum image obtained from the Hi data in 
March 2002 is shown in Fig. |6l This and the following im- 
ages have been rotated by 4? 9 counter-clockwise in order to 
optimise the information in optical depth slices. 

W e retrieved the data presented by iGreisen fc Liszt I 
1 1986 ) from the VLA archive and reduced them in order 
to check whether our results were consistent. The total con- 
tinuum flux density detected by MERLIN in the 10-arcsec'^ 
region around component '1', shown in Fig. (6] is 10.5 Jy, 
compared to a peak flux density of 12.1 Jy in the newly 
constructed VLA image at (5.7 x 4.2) arcsec'^ resolution. 
Thus, the MERLIN continuum image contains almost 90 per 
cent of the flux density as observed at lower resolution over 
the relevant field of view. The expectation that the MERLIN 
line data represents most of the VLA emission/absorption in 
this region is verified by the comparison of the two integrated 
H I spectra over the extent of the continuum source shown in 
Fig. |6] The VLA and MERLIN data agree remarkably well, 
given the uncertainties. In addition, the MERLIN spectrum 
over a much smaller region of (0.35 x 0.35) arcsec^, shown in 
FiglZl fl'lso agrees with the mean VLA spectrum as expected. 
There remains the possibility that for 3C 161 the opacities at 
T > 1.3 may be somewhat overestimated (by 5t ~ 0.2 — 0.3), 
if there is a slight amount of missing continuum flux density. 
The weak extended components of 3C 161 (see Fig |A4|l are 
not detected in the MERLIN data due to the larger angular 
sizes and resulting lower surface brightnesses (giving < 0.1 
mjy beam~^ on the scale of the MERLIN beam) compared 
with the bright NE lobe. There is some evidence in Fig. [7] 
for a separate line at « 7 km s~^, close to a line de tected in 
the ea rlier NRAO 3-element interferometer data of [GreisenJ 
i|l976t) 

The observed range of velocities from to 35 km s^^ 
suggests the infiuence of Galactic rotation in the third quad- 
rant of the Galaxy with more positive velocities arising from 
IS gas at larger distances from the Sun. As an example the 
velocity of 34.2 km s~^ would have an inferred distance of 
~ 3.5 kpc. The distances that can be associated with gas at 
lower velocities ( less than ~ 12 km s~^, twice the cloud-to- 
cloud velocity dispersion) are quite uncertain and a better 
estimate may well be to use Equation [1] w ith a resultant 
distance of 0.8 kpc. (jGreisen fc Liszt I (|l986l ) have discussed 
the problem of distances of Hi features in this direction. 
They point out that for 3C161 : " ... the Gould Belt phe- 
nomenon carries the cloud velocities to rotationally permit- 
ted values.." ) 

The MERLIN beam has an NS:EW axial ratio of ~ 
3.5 : 1 and the resolved emission from 3C 161 is effectively 
only one beamwidth wide in the N-S direction. We therefore 
show slices in the direction of the long axis of the source, in 
the direction of the most favourable resolution. Figures [S] 
[HI 1101 [TT]and [T2]show optical depth images with inserts 
showing slices of H I optical depth variation in single chan- 
nels at velocities of 4.5, 6.6, 11.5, 20.2 and 34.2 km s~^. In 
the analysis following, we will derive distances for the lower 
velocity H I features (Figs. [H] [9] and I10[) by assuming that 
they arise from gas at 0.8 kpc, where the MERLIN beam of 
0.13 arcsec corresponds to about 100 AU. We will use the 
flat Galactic rotation curve for the higher velocities. These 
estimates are all quite uncertain. 



The variation of t(H i) at 4.5 km s~^ (Fig. [8]) is complex. 
The strongest H I absorption, at a positional offset along the 
slice of 1.8 arcsec, shows an excess opacity (with respect to 
the local background opacity) at a level of ~ 1. This feature 
is spatially resolved with an angular size of ~ 0.5 arcsec, or 
400 AU. Its signal-to-noise ratio is > 30a". The other fea- 
tures in the slice show opacity fluctuations with poor signal 
to noise which are probably not significant. At a velocity 
of 6.6 km s~^ (Fig. ^, there is a prominent Hi peak at a 
displacement of 1.3 arc sec; t(Hi) fluctuates by ~ 1.4 over 
a region of angular size < 0.15 arcsec (120 AU). These two 
spectral features overlap but can be decomposed into Gaus- 
sian hue profiles, each of FWHM 2.3±0.2kms"\ 

Farther East, the Hi distribution is almost constant at 
r(Hl)=0.8 over an angular extent of at least 0.7 arcsec. 
Fig. \TU\ shows two prominent peaks at a velocity of 11.5 
km s~^ at angular displacements of 1.8 and 1.4 arcsec, with 
excesses in t(Hi) of 0.5 and 0.6 and angular sizes of 0.5 and 
0.2 arcsec (400 and 160 AU), respectively. The correspond- 
ing spectral feature has a FWHM of 3.2 ± 0.3 km s^^. 

The Hi opacity at 20.2 km s'^ (Fig. [TT]) is an order 
of magnitude weaker than at lower velocities. This spectral 
feature has a rather uneven shape with a FWHM of 2.5 ±0.5 
km s~^. The slice shows at least one Hi component at an 
offset of 1.7 arcsec with fluctuations in ■7"(Hl) of ~ 0.1 and 
an angular size of 0.5 arcsec (900 AU at a distance of 1.8 
kpc). The highest velocity and hence most distant Hi, at 
34.2 km s~^ fFig. I12|) . has a well-defined spectral profile of 
FWHM 3.0±0.2 km s"\ The line does not show evidence of 
any significant variation in the Hi opacity over an angular 
scale of ~ 1 arcsec (3500 AU at a distance of 3.5 kpc). 



4 SUMMARY AND CONCLUSIONS 

We have used the MERLIN array to image Galactic H I ab- 
sorption, with a resolution of ~ 0.1 arcsec, in the direction 
of three background sources, 3C 111, 3C 123 and 3C 161. 
We resolve the Hi distribution towards the latter two ob- 
jects with angular sizes ^ 1 arcsec. In the case of 3C 161, 
prominent spatial variations at velocities close to zero (in 
the range 4.5-11.5 km s~^) have been observed. In addi- 
tion, at 20.2 km s~^ the opacity of the Hi is an order of 
magnitude lower (« 0.1 — 0.2); at this velocity significant 
fluctuations in the Hi opacity across the face of 3C161 are 
observed. The inferred distance of this Hi feature is 1.8 kpc. 
At the highest velocity in this direction (34.2 km s~^ at an 
inferred distance of 3.5 kpc) no significant changes in the 
weak Hi features with r(Hl)~ 0.1 are observed. In the di- 
rection of 3C123, there are suggestions of opacity variations 
in the low-velocity Hi, but with a much lower statistical 

significance. 

As pointed out bv lGreisen fc Liszt I l|l986r ), the expec- 
tation would be that the Hi features at higher velocities 
would subtend smaller solid angles and thus cover a smaller 
fraction of the continuum background source; however the 
existence of a cloud size spectrum and random cloud veloci- 
ties could mask the expected behaviour in this small sample. 
In addition, the higher- velocity Hi lines have smaller opac- 
ities (0.1 to 0.2) compared to the lower velocity lines (r ~ 
1). This decrease in mean opacity reduces the sensitivity 
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of these studies with respect to detecting significant spatial 
variations in the higher velocity lines. 

On a spatial scale of 400 to 500 AU, prominent changes 
in opacity a t the level of ^r ~ 1 in the direction of 3C 161 
are striking. iBrogan et al.l (|2005l ) had previously determined 
5t ^ 0.5 with a linear resolution of ~ 25 AU in the direc- 
tion of 3C 138. The peak change in 5t of about 1 seen in 
the MERLIN data is consistent with the peak fl uctuation of 
about 0.7 on scales of 400-500 AU predicted bv lDeshpandel 
l|200CI ). Previous studies of small scale Hi have tradition- 
ally calculated the implied densities based on the assump- 
tion that the line of sight dimension is comparable to the 
observed transverse dimension. iHeilesI ((20071) has discussed 
this assumption in detail; the expectation is that the ob- 
jects are far from isotropic. The 400-500 AU structures de- 
tected in these observations would imply volume densities 
of < 5 X 10^ cm~^ with column densities a f ew times 10^° 
cm~^ if the excitation temperature is 50 K IjBrogan et al.l 
120071 ) and the objects have a path length along the line of 
sight that is more than ten times the transverse dimension. 
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Table 1. Source and Observing Parameters. The duration is the total length of each observation, about 2/3 of which was spent on the 
main target. The RMS noise crrms given for the line observations is for a single spectral channel. 



Source 




3C 111 


3C 123 


3C 161 


Galactic Coordinates 


(1°, b°) 


161.68, -8.82 


170.58, -11.66 


215.44, -8.07 


Broad-band Continuum 










Date 


YYYMMDD 


20030202 


20030204 


20011030 


Duration 


(hr) 


17.5 


16 


10 


Frequency 


(MHz) 


1408 


1408 


4994 


Phase Reference Source 




B0415+398 


B0430+289 


B0605-085 


Target Peak 


(Jy bm-i) 


1.285 


0.749 


0.539 


RMS noise 


(mjy bm^-') 


2.0 


2.6 


0.6 


Narrow-band/Line 










Date 


YYYMMDD 


20030508 


20030501 


20020304 


Duration 


(hr) 


16 


16 


9 


Vlsr 


(km s-l) 


5.0 


3.0 


17.5 


Phase Reference Source 




B0415+398 


B0436+306 


B0605-085 


Phase Ref. Flux Density 


(Jy bm-l) 


0.41 


0.55 


2.1 


Synthesized Beam 


(arcsccxarcsec) 


0.19 X 0.15 


0.23 X 0.14 


0.46 X 0.13 


Position Angle 


(deg) 


45 


28 


21 


Continuum Peak 


(Jy bm-l) 


1.27 


0.75 


3.40 


Continuum crmis 


(mJy bm~^) 


2 


4 


4 


Line CTrms 


(mJy bm^^) 


8.5 


11 


18 
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Figure 1. The Hi optical depth spectra towards the core and the resolved NE component of 3C111 are shown by the thin black line 
and the thick grey line, respectively. The velocity channel separation is 0.41 km s~^ and velocities are given with respect to the Local 
Standard of Rest. The optical depth was averaged over areas of 0.3 arcsec X 0.3 arcsec for the core and 0.45 arcscc X 0.65 arcsec for the 
NE component. 
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Figure 2. MERLIN continuum image of 3C 123 observed at 1408 MHz, using 13.5-MHz bandwidth. The beam is 0.23 arcsec X 0.14 
arcsec at a position angle of 28°, indicated in the lower right corner. The contour levels are (-3, 3, 6, 12, 24, 48, 96, 192) X the rms noise 
level of 2.6 mjy beam^^. The weak core is detected as well as the two prominent hotspots. The greyscale shows the flux density in mjy 
beam^^. 
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Figure 3. Mean Hi optical depth spectrum towards the peak emission component (SE) in 3C 123 as observed with MERLIN, averaged 
over an area of 0.25 arcsec X 0.25 arcsec centred at 041^ 37™ 05?01 +29°40'11''0 (J 2000). The velocity charmel separation is 0.41 km 
s~^ and velocities are given with respect to the Local Standard of Rest. 



10 Goss et al. 


29 40 12.5 
12.0 
11.5 
11.0 
10.5 
10.0 
09.5 
09.0 
08.5 
08.0 



o 
o 
o 



o 

c 
"o 

0) 
Q 




0.5 1 1.5 

Dec. offset from 29° 40' 09". 57 (csec) 



"T^ 



1.5 1 0.5 

RA offset from 04^ 37'^ 04^917 (osec) ■ 






04 37 05.1 



05.0 04.9 04.8 04.7 

Right Ascension (J2000) 



04.6 



04.5 



Figure 4. MERLIN Hi optical depth image towards tlic SE component of 3C 123 in a single velocity channel centred at 4.2 km s~ , 
adopted distance to Hi 530 pc. The greyscale represents Hi optical depth, and the contours represent continuum emission at (-3, 3, 6, 
12, 24, 48, 96) x the rms noise level of 6 mjy beam~^ for 0.15-MHz bandwidth of line- free channels at 1420 MHz. The Hi optical depth 
was not calculated where the continuum is < 3 per cent of the peak or where the H I line is completely saturated. The inserts show slices 
in Hi optical depth along the white portions of the horizontal and vertical lines, from E-W (West at right, lower plot) and N-S (south 
at left, upper plot) respectively. Symmetric ilcr error bars are indicated; note that at the highest opacities the upper uncertainties may 
be greater than the lower. 
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Figure 5. Hi optical deptli image towards 3C 123 in a single channel centred at -20.1 km s ^, adopted distance to Hi >530 pc. The 
insert shows an E-W slice. See Fig. |4] for more details. 
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Figure 6. MERLIN continuum image of 3C 161 observed at 1420 MHz in 0.18-MHz bandwidth of line-free channels. The image has 
been rotated by 4?9 counter-clockwise. The contour levels are (—3, 3, 6, 12, 24, 48, 96, 192, 384) X the rms noise level of 4 mjy beam"'^. 
The beam is 0.46 arcsec X 0.13 arcsec at a position angle of 21°, indicated in the lower right corner. The greyscale shows the flux density 
in mJy beam^^. 
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Figure 7. The Hi optical depth spectrum towards the continuum peak of 3C 161, averaged over 0.35 arcsec X 0.35 arcsec. The velocity 
channel separation is 0.41 km s^^ and velocities are given with respect to the Local Standard of Rest. Note the wide range of velocities 
over which the H I absorption is observed. 
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Figure 8. MERLIN Hi optical depth image toward 3C 161 in a single channel centred at 4.5 km s~^, adopted distance to Hi 0.8 kpc. 
The greyscale represents H I optical depth, and the contours represent continuum emission as in Fig. \E\ The H I optical depth was not 
calculated where the continuum is < 1 per cent of the peak. The image has been rotated by 4? 9 counterclockwise. The insert shows a 
slice in Hi optical depth along the white portion of the horizontal line (West at right). Symmetric ±1(t error bars are indicated; in some 
cases these are less than the width of the line. Note that at the highest opacities (r > 1.5), the upper uncertainty is likely to be 1.5 to 2 
times larger than the lower uncertainty (which is plotted). 
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Figure 9. MERLIN Hi optical depth image toward 3C 161 in a single channel centred at 6.6 km s~^ , adopted distance to Hi 0.8 kpc. 
See Fig.Ofor more details. The opacities shown in the insert for the compact component at an offset of 1.3 arcsec may be overestimated 
(by 0.2 - 0.3), if there is a slight amount of missing flux density in the continuum. 
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Figure 10. MERLIN Hi optical depth image toward 3C 161 in a single channel centred at 11.5 km s ^, adopted distance to Hi O.i 
kpc. See Fig. [8] for more details. 
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Figure 11. MERLIN Hi optical depth image toward 3C 161 in a single channel centred at 20.2 km s ^, adopted distance to Hi l.i 
kpc. See Fig. [8] for more details. 
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Figure 12. MERLIN Hi optical depth image toward 3C 161 in a single channel centred at 34.2 km s ^, adopted distance to Hi 3.5 
kpc. See Fig. [8] for more details. 
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APPENDIX A: 
OF 3C161 



CONTINUUM OBSERVATIONS 



We describe the radio continuum images of 3C161 obtained 
as part of this work since they are not published fully else- 
where. Although a powerful radio source, 3C161 lies at low 
Galactic latitude (6 = —8.07°) and has therefore escaped 
much attention si nce it has often been ex cluded from sam- 
ples under study. iGreisen fc Liszt I l|l986t ) showed that this 
source possesses a classical double radio structure with an 
angular size of 62 arcseconds, albeit with a large brightness 
ratio in excess of 60 between the Northern and Southern 
hotspots. The H I absorption reported in this paper is asso- 
ciated with the brighter Northern hotspot. 

A MERLIN high-resolution 5-GHz image is shown in 
Fig. lAll The observations are summarised in Table [l] and 
the data reduction is outlined in Section[2l The point source 
calibrator, B0552-I-398, was found to have a flux density of 
5.74 Jy by comparison with 3C286. 3C 161 and the nearby 
(5.5°) phase-reference source B0605-085 were observed alter- 
nately for 5.5 and 2 min, respectively, over a total of 9.5 hr. 
We also used B0605-085 to correct for polarization leakage 
and we calibrated the polarization angle (P.A.) with respect 
to 3C286 which has an assumed P.A. of -1-33° (in the map 
plane) . 

Figure IXT1 shows a conventional hotspot structure (peak 
539 mjy beam^^). The total intensity map noise around 
the hotspot is increased about three-fold with respect to 
more distant regions of the map (and the polarized intensity 
images) due to sidelobes of this bright source. 

The polarized intensity image, 'Fm. \A2\ { left], shows that 
the brightest polarized emission lies around the outer edge 
of the hotspot region. The maximum polarized intensity de- 
tected is 145 mJy beam~^. The plotted vectors indicate the 
directions of the in situ E field, which have been corrected 
for foreground rotation us ing the rotation measur e of 112 
rad m^ given bv Vallee. Broten fc MacLeodI l|l988l ). These 
vectors indicate a circumferential magnetic field, as would 
be expected around the outer edge of the working surface of 
a lobe. The trajectory of the incoming beam or jet can be 
traced in the polarized intensity image. Fig. IA2I (left). The 
jet enters the hotspot region from the West where the E 
vectors in Fig. lAll are parallel to each other and perpendic- 
ular to the jet direction (implying a magnetic field aligned 
with the jet) . The fractional polarization is highest here and 
at the leading edge of the lobe, at ~ 30 per cent. The po- 
larized intensity falls to nearly zero just as the jet enters 
the hotspot in a region where the polarized position angle 
is rotating rapidly, presumably due to shocks. 

In addition to the bright emission associated with the 
Northern hotspot, we detect a compact component of 1.15 ± 
0.12 mJy beam~^ (lOcr at the local rms noise level) at a 
distance of 25.5 arcsec SW of the hotspot, shown in Fig. lA2l 
(right). Its position is J 2000 right ascension 06''27™09'189, 
declination — 05°53'26'.'43 (uncertainty 25 mas). This is a 
candidate for the core of 3C 161; however, we note that it is 
not seen in any other available data. MERLIN did not detect 
any emission from the much weaker Southern hotspot. 

Figure IA3I shows an image of 3C 161 at 5 GHz which 
we made from VLA archive data (proposal code AP370, 
observed on 19980915 in the B-array configuration). The 
source was observed as a calibrator and, to our knowledge. 



no maps at this frequency and resolution range have previ- 
ously been published. The positions of the MERLIN hotspot 
and core candidate are marked. The VLA did not detect any 
core emission but the 3a noise level of this image is close to 
the peak emission of the core candidate measured by MER- 
LIN, so it is feasible that it was just below the detection 
threshold of these VLA data. 

Figure IA4I shows the 1.4 -GHz continuum count erpart 
to the Hi data published bv iGreisen fc Liszt] l|l98d ). The 
data have been re-processed and the image shown here is 
naturally-weighted. This map is thus somewhat deeper and 
more sensitive to extended lo w-surface-brightn e ss em ission 
than the original, figure 2a in IGreisen fc Liszt I (|l986l ). The 
Northern component of this FRII double source is a little 
unusual since it has a tail which arcs around through ~ 
180 degrees to the SE. This is in the same sense as the jet 
entry direction into the hotspot. There is also low surface 
brightness emission SW of the double structure which may 
be associated with the source. 

A linear feature, seen in Fig. IA4I extending from NW 
of the Northern hotspot, is difficult to interpret. It could be 
associated with a faint optical object seen in the Digitised 
Sky Survejo which also lies just NW of the 3C161 hotspot, 
but this cannot be related to 3C161. Indeed, there are no 
obvious optical counterparts which can be related to any 
part of the observed radio structure of 3C161, and thus it 
remains unidentified. 



^ |http : //archive . stsci . edu/dss/index .html | 
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Figure Al. MERLIN continuum image of 3C 161 observed at 4994 MHz using 15-MHz bandwidth. The contour levels are (—3, 3, 6, 
12, 24 ...) X the rms noise level of 0.35 mjy beam~^. The circular beam of 100 mas FWHM is indicated in the lower right corner. The 
greyscale shows the polarized intensity in mJy beani^^. The vectors represent the direction of the E field and 10 mas length corresponds 
to 4 mJy beam^^ polarized flux density. 
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Figure A2. MERLIN continuum images of 3C 161 observed at 4994 MHz using 15-MHz bandwidth, left Polarized intensity. The contour 
levels are (-1,1, 2, 4 ... 32) X 4 mjy beam~^. The circular beam of 100 mas FWHM is indicated in the lower right corner, right The core 
candidate in total intensity. The contour levels are (-3, 3, 6, 9) X the local rms noise level of 0.12 mJy beam^^. The beam is 103 mas 
x47 mas at a position angle of 19°, indicated in the lower right corner. 
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Figure A3. VLA continuum image of 3C 161 observed at 4860 MHz using 100-MHz bandwidth. The contour levels are (—3, 3, 6, 12, 24 
...) X the rms noise level of 0.33 mjy beam~^. The beam is 4'.' 94 x 2'.'71 at a position angle of -36°, indicated in the lower right corner. 
The crosses mark the positions of the MERLIN hotspot peak and the core candidate. 
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Figure A4. VLA continuum image of 3C 161 observed at 1420 MHz using 0.15-MHz bandwidth. The contour levels are (-3, 3, 6, 12, 
24 ... ) X the rms noise level of 2 mjy beam"'^. The beam is 7'.'19 x S'.'Ol at a position angle of 17°, indicated in the lower right corner. 



